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Abstract 

It is necessary for the statistical description of collective effects in liquids to 
set that or other approximation between direct and pair correlation functions 
which describe a neighboring order. The analytical solution of the generalized 
Ornstein-Zernike (OZ) equation was obtained for a system of particles at 
T — 0. It is shown that in this limit case the neighboring order disappears 
and the direct correlation function describes a distant order which is typical for 
the ideal crystal. The approximation correctly describing the limit transition 
liquid-solid at T = will have a physical meaning. 

Introduction 

The structure and thermodynamics of both dense gases and liquids 
is described by a pair correlation function hu considering the collective 
effects in dense mediums. Nowadays about 20 approximate integral 
equations of the Ornstein-Zernike (OZ) type connecting the direct Cm 
and pair hu correlation functions are proposed for describing of the 
above-mentioned one |T], g]. The degree of exactness of every such 
approximation is impossible to be evaluated. The physical meaning of 
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their basic approximations has not been cleared up to the end yet. Thus 
such a situation suggests a choice of maximum exact approximations 
in a comparison with the result achieved at the numerical experiment. 
Nevertheless such a comparison doesn't reveal the physical meaning of 
approximations. 

The OZ equation is to be used for the spatial homogenous and 
isotropic mediums as in such a case a two-particle distribution func- 
tion (as well as a pair correlation function) depends only on a mutual 
distance between the particle centers 7*12, that makes all the calculations 
much more simplified. 

But there exists a more generalized form of the OZ equation for 
the spatial inhomogenious anisotropic mediums 0, §], [5|]. It has not 
been used widely yet for the solution of more exact problems. One 
can mention its use only for the description of the liquid-solid phase 
transition [5|, |6j] . This work describes the use of the Martynov-Sarkisov 
approximation, which gives the best agreement with the data of the 
numerical experiment for the hard spheres system |5j. 

In this work we suggest using the generalized OZ equation for the de- 
scription of the ideal crystal at T = 0. The pair correlation function h\2 
in this limit case takes the meaning of the Dirak 8 function. As a result, 
a linear integral equation for the direct correlation function is achieved, 
which has a simple analytical solution. It should fit The results known 
from crystallophysics for the ideal crystal. So the limit transition to the 
model of the ideal crystal may be considered as a physical criterion for 
the evaluation of the exactness for the approximations used in physics 
of liquids. 

1 Initial equations 

In works [|3|, [5], [7j it was shown that a hierarchy of Bogolubov- Born- 
Green-Kirkwood-Yvon (BBGKY) equations may be modified into a two 
equations system 




(i) 




(2) 



2 



for two unknown functions: a one-particle distribution function G\{fi) = 
exp(co>i(fi)) and a two-particle distribution function 

Gijifi^^GiG^l + hifr (3) 

hij = -1 + exp(-^ + Slij(fi, fj)), (4) 

where n = N/V is density, -molecular interaction potential, = kT- 
temperature, d(i) = dfi- a volume integral element of the i particle. 

All previous distribution functions are expressed through G\ and Gyi- 
The direct correlation functions Sij, Cij in equations (p and (@) are 
evidently expressed through infinite integral sets distinguished from the 
products of the pair correlation functions 

S tJ = h i:j - Mj - X -h l3 (% + i Afg> ) , (5) 

C ij = h ij -n ij + ^M§\ (6) 

in which 

M$ = n 2 J J G z Gihahuhuh 3j h4jd(3)d(4:) + k = 1,2 (7) 

the so-called bridge-functionals represent infinite sums of non-reduced 
diagrams. 

For spatial homogeneous and isotropic systems (dense gases and liq- 
uids) at u\ = 0, G\ — l,Gi2 = ^12(^12), equation ([I]) is a definition of 
the activity factor logarithm In a. Equation (g) describing the neigh- 
boring order in liquids tends to come to a well-known OZ ratio 

h 12 = C 12 + nJC 13 h 13 d{3) : (8) 

being the basis of the modern theory of liquids and molds |2j . For 
making ratio (|8|) closed it is necessary to set an approximation connect- 
ing the direct correlation function C^iif 12) with the pair correlation 
function ^12(^12). 
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2 The ideal crystal 



One particle distribution function has the meaning of local density 
in laboratory coordinate system. For spatially homogeneous isotropic 
mediums one particle distribution function equals one because the ex- 
ternal field does not exist. For the crystal condition one particle distri- 
bution function does not equal one, even if the external field does not 
exist. To describe such a system it is necessary to introduce the ex- 
ternal field according to Bogolubov which fixes the crystal place in the 
space and then to make it tend to zero. Then one particle distribution 
function describes the local density distribution in the regard to the 
fixed coordinate system giving the crystal place as a whole. As an ori- 
gin of such a coordinate system it is convenient to choose one of lattice 
particles giving it number zero conditionally. Therefore one particle 
distribution function G\ = exp(o; 1 ) for the crystal can be represented 
in the form of the periodic function 

Gi(f 1 ) = 'E%exp^r 1 ), (9) 

k 

For the ideal crystal at T = every Fourier's component G% is similar 
and G\ is a superposition of three-dimension Dirak 5 functions 

Gi(ri) = E*(ri-fy (10 ) 

where the summation is done by crystal lattice nodes f n . whose coordi- 
nates, in their turn, must be defined. 

The mutual arrangement of lattice nodes (that is the lattice period 
and the type of the crystal system) is defined by the potential of molec- 
ular interaction. On the other hand, the potential $12 enters the pair 
correlation function in the form of the combination — Let us con- 
sider the kind of the function hu at T = 0. Let the molecular interaction 
be described by Lenard- Jones potential 

*i2 = 4 £ [(-) 12 -(-f], (11) 

where a is a characteristic size of a molecule. The potential 3>i2 has a 
minimum at the point tq = a^2. Therefor the function — ^ at T — > 
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has S -shaped maximum at this point. That is why one can neglect the 
thermal potential U12 describing the neighbouring order. As a result 
hu is approximated one-dimension Dirak function 5(ru — Vq). 

As a result at T = we get the following equation for the straight 
correlation line with regard to (|TU1) 

C[f +nEC 1 ( 2 ) (f 1 ,f n )^(r 2n - r ) = S(r l2 - r„). (12) 

We notice that (jT2|) is the linear equation whose solution is sought in 
the form 

Cu = a5(r u - r ), (13) 

where a is a coefficient to be defined. Substitution C\ 2 in (|T2|) results 
in the following expression 

S(rn-r ) = a£(ri 2 -r )+na^£(| ?i-r n \ -r )6(\ r 2 -r n | -r ). (14) 

Particles 1 and 2 are the neighbouring ones. The summation is done 
by those nodes which are the nearest both for the first particle and for 
the second one. Let such particles have No numbers. Using fll^p we get 

ot = - ■ 1 , r , r u = r , ri3 = r , r 23 = r (15) 
1 + nA 

The correlation defines the crystal structure (the lattice period and 
the type of the crystal system). For determining the crystal system we 
should calculate the number Nq. 

At T = the system is densely packed. There are two types of 
lattices with the dense packing - face centered cubic and hexagonal 
lattices [|J. Their structure is formed by successive layers of densety 
packed planes. Every particle has twelve neighbouring particles. This 
structure type is connected with the distant order in the crystal. But 
in any substance (including the crystal) there is the neighbouring order 
which is described by the second addend in the right-hand side (|T|) in 
this case. 

Let us take arbitrarily a couple of the neighbourong particles ar- 
ranged in the plane (111). It has four particles juxtaposed with it for 
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both types of the lattice: two particles in the plane (111) and one parti- 
cle in every adjoined plane. If particles 1 and 2 lie in two neighbouring 
planes the situation is different. For the face centered cubic lattice Nq 
is 4 as before but for the hexagonal lattice NQ is 3. So Nq =const for 
the face centered cubic lattice. As the pair of the neighbouring particles 
is chosen arbitrarily this indicates that the face centered cubic lattice 
is realized. It corresponds to the experimental data for crystals of inert 
gases. However for Lenard- Jones potential the hexagonal structure has 
much lower energy than the face centered cubic lattice therefor it is this 
lattice that must be realized. The problem of the crystal system is not 
solved in this paper yet. 

Now we show that solutions ( |IUD and (|I3|) satisfy the equation ([ID. 
The direct correlation function neglecting two particle correlations 
UJ12 equals 

C[l ] = c{f = a5(r 12 -r ). (16) 
Substitution this expression in ([T]) we shall get 

uji = na J2S(\ n - f n | -r ) + In a. (17) 

Let a = 1 and doing summation by fu as it is done in ( \PSj ) 

uj\ = naN()5(\ r\ — Rq \ —vq) + In a. (18) 

— > 

where Rq is a radius of the vector-particle placed in the coordinate 
system origin. Since coordinates of all the nodes are known now, u\ is 
a given function G\ can be given in the form (|I0|) . 

The complicity of problems being solved by the physics of condensed 
matter makes one use different physically based approximations. In 
particular, in the physics of plasma and solids an approximation of the 
self-coordinated field is widely used. 

In the physics of liquids a superpositional approximation is well- 
known, which is used to express a triple function of the distribution 
through two-particle ones. It is to the superpositional approximation, 
that fits to the hypernetted chain equation, which results from (@) by 
the neglect of all the irreducible diagrams. Summing up of such dia- 
grams is practically unrealizable. Physically it is connected with the 
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fact that at the limit temperatures it is necessary to take into account 
the collective effects, stimulated by a molecular motion. That is why 
it is necessary to offer a physical value criterion of the irreducible dia- 
grams. 

A limit transition to the ideal crystal at T = represent such a 
criterion. In this case it is possible to neglect the collective effects 
connected with the molecular motion. As a result a limit ratio for the 
direct correlation function Cu, that effectively includes a contribution 
of the irreducible diagrams, is possible to be obtained. 

From a great number approximations between Cn and hyi the most 
physically valid is that one which fits to the given limit ratio. 

Conclusion 

The solution of Ornstein-Zernike generalized equation is obtained 
for the ideal crystal at T = 0. In this case the pair correlation func- 
tion h\2 Degenerates to Dirak S function and for the direct correlation 
function Cyx the linear integral equation is found which has the ana- 
lytical solution. By using it the main problem of physics of liquids is 
solved - establishment of closure between the direct and pair correlation 
functions. 

At the same time there is still a problem for the aggregate condition, 
the method developed by us can be used as a criterion for evaluating the 
closures: those resulting in the solution for the ideal crystal at T = 
will have the physical meaning. 
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